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Abstract 

This paper presents synthesis methods for the de- 
sign of constant-gain positive real controllers for pas- 
sive systems. The results presented in this paper, in 
conjunction with the previous work by the authors on 
passification of non-passive systems, offer a useful syn- 
thesis tool for the design of passivity-based robust con- 
trollers for non-passive systems as well. Two synthesis 
approaches are given for minimizing an LQ-type per- 
formance index, resulting in optimal controller gains. 
Two separate algorithms, one for each of these ap- 
proaches, are given. The synthesis techniques are 
demonstrated using two numerical examples: control 
of a flexible structure and longitudinal control of a 
fighter aircraft. 

Introduction 

Passivity-based controllers have been proved to be 
highly effective in the control of inherently passive sys- 
tems [Jos. 89, Kel.96]. Recently, it has been shown that 
these controllers can also be used for control of non- 
passive systems as well once such systems are passi- 
fied by techniques introduced in [Kel.97, Kel.98]. The 
main advantage of using such controllers is the sta- 
bility robustness of the closed-loop system. There 
are numerous results available in the literature on the 
passivity-based controllers. However, most of these re- 
sults have focused on the analysis part and not much 
work has been done on the synthesis of such con- 
trollers. The limited results available to date on the 
synthesis of positive-real controllers can be found in 
[Loz.90, Had. 94, Saf.87]. In [Loz.90] an LQG-based 
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design technique was given whereas in [Had. 94] an 
7f 2 /Roo-based design procedure was given. Both of 
these methods used stochastic models and are compu- 
tationally intensive. 

This paper gives two approaches to synthesize 
constant-gain positive-real (PR) controllers. The first 
approach is based on the use of symmetric gain matrix 
whereas the second approach allows a non-symmetric 
gain matrix. In the case of the first approach, the 
necessary conditions are derived by minimization of a 
suitable LQ performance index. A synthesis algorithm 
is derived based on these necessary conditions. For the 
second approach, a modified version of the algorithm 
proposed in [Moe.85] is used. The organization of the 
paper is as follows. First, we present the problem for- 
mulation followed by two separate algorithms to com- 
pute optimal controller gains, and finally, two numer- 
ical examples to demonstrate the synthesis methods. 

Positive-Real Optimal Controllers 

This section gives the formulation of optimal con- 
trol problem wherein it is desired to synthesize a 
constant- gain strictly PR (SPR) output feedback con- 
troller for a Linear Time Invariant (LTI) PR system. 
Consider a Positive-Real (PR) LTI system: 

x — Ax + Bu (1) 

y-Cx (2) 

where x(t) E 3£ n , u(t) E y(t) E K m , and input 
u(t) is given by the output-feedback control law: 

u — —Gy (G = G t > 0). (3) 

The plant (1) being PR satisfies the following con- 

straints as a result of Kalman- Yakubovich lemma for 


1 



some positive definite P = P T G 5J nxn , and L E 

r xn . 

A t P + PA = -L t L (4) 

C = B T P (5) 


Assuming x(0) to be a random vector, J in Eq. (13) 
can be replaced by its expected value, and thus Xo in 
Eq. (14) denotes £[z(0)£ T (0)]. 

The optimization problem stated above can now 
be re- phrased as follows: 


Let us suppose that it is required to find an asymptot- 
ically stabilizing constant- gain, strictly positive-real 
output feedback controller G for the closed-loop sys- 
tem given by Eqs. (l)-(3) such that the following per- 
formance index is minimized. 

/»oo 

J == / ( x t Qx + u T Ru)dt (6) 

Jo 

This is an optimization problem where the perfor- 
mance index (6) is to be minimized subject to the 
constraints (l)-(5). Using control law of Eq. (3) the 
closed-loop system becomes 

x = {A-BGC)x. (7) 

The constraint that G is SPR can be imposed by forc- 
ing G to satisfy 

G = KK t (8) 

where K E 3? mxm has full rank. Substituting Eq. (8) 
in Eq. (7) yields 

x — (A - BKK t C)x = A d x (9) 

where A c i — A — BKK T C . The performance function 
(6) can be rewritten using Eq. (3) as 
1*00 

J= ( x T (Q + C T K T I<RKK T C)x)dt . (10) 

Jo 

Further, it can be shown that 

J= ^ T (0)£x(0). (I 1 ) 

where E satisfies: 

A T d E + E A ci + Q + C t K t KRKI< t C = 0. (12) 

For a given K, R> and Q this is a Lyapunov equa- 
tion to be solved for E. Now using the trace iden- 
tity, tr(AB)=tr(iM), the performance index J, can 
be rewritten as 

J = tr(EXo) (13) 

where 

(14) 


Minimize 
Ky E 

subject to 


J = tr(EXo) (15) 

g ~ A d T Yl + E A d -f- Q — 0 (16) 


where Q = Q + C T G T RGC . The problem defined 
by Eqs. (15-16) is a constrained minimization problem 
which can be further modified into an unconstrained 
minimization problem by using Lagrange multipliers 
and augmented performance function. 

Let S = S T be the Lagrange multiplier with con- 
sistent dimensions. If Hamiltonian is defined as: 


H = tr(EA 0 + gS) = tr(EAo) 

+ tr ([A - BKK t C] t 1;S + E [A - BKI< T C]S 

+ QS + C t KK t RKK t CS^J (17) 

the constrained optimum of Eq. (15) subject to con- 
straint (16) is same as the unconstrained optimum of 
Eq. (17). The necessary conditions for an optimum 
are then given by setting the partial derivatives of H 
with respect to its arguments to zero, i. e., 

=> g = ^4 c / T E -j- E A d + Q = 0 (18) 

=> Jrrtr ( - [C' t KK t B t T,S + E BKK T CS] 

+ C t KI< t RKK t CS ) = 0 

=► - [CS T e t b + b t e sc t + b t e t s t c t 

+ CST,B ] K + 2[CSC t KK t R 
+ RKI< t CSC t ]K = 0 
=> [CSC t KK t R + RKK t CSC t - CSEB 
- B t HSC T ]K — 0 (19) 

^ Xq + A d S + SA cl T = 0 (20) 

Note that the equation (18) is just the constraint 
equation, whereas Eqs. (19) and (20) give necessary 
conditions for an optimum. The unconstrained opti- 
mum is a local minimum if the Hessian H KK is posi- 
tive definite. 


^ = 0 
dS 

d H _ 

dK 


0 


dH 

as 


Xq = x(0)a.’ T (0). 



The following theorem gives conditions under 
which the the performance function (15) decreases, 

and leads to a synthesis algorithm. 

Theorem 1- If E 1 , E°, K 1 , and K° are solutions to 
the following equations 

A c or E° + £°A C ° + Q° = 0 (21) 

Ac^S 1 + £ 1 A C 1 +Q 1 = 0 (22) 

Xo + A C °S° + S°A c 0T = 0 (23) 

where: A i c - A- BG’C = A-BK'K^C and Q* - 
Qi + C T G iT RG'C = Q i + C T K i K iT RK' K iT C and 
if K 1 satisfies 

CS°C T K 1 K 1T R + RK 1 K 1T CS°C T - 

CS°'E 1 B - B t 'S 1 S 0 C t = 0 (24) 

then 

J{K°) — J(K l ) > 0. (25) 

provided that A l c is Hurwitz for i = 0, L 

Proof- The proof is omitted due to space limitations. 

Using the necessary conditions (Eqs. 19 and 20) 
and Theorem 1, a numerical algorithm is presented 
below which can be used to synthesize the PR gain 
matrix G. 

Synthesis of symmetric PR gain: 
Algorithm 1 

Using Theorem 1, the following iterative algorithm 
can be obtained to compute the gain G. 

Step 1. Choose A" 0 and solve for D° using equation 

( 21 ). 

Step 2. Solve for S° using equation (23). 

Step 3. Solve for K l and X 1 simutanously by using 
equations (22) and (24). 

Step 4. Set A" 0 = A 1 , go back to step 1. 

Iterate until convergence is obtained. 

It should be noted that Step 3, which consists of 
solving nonlinear coupled matrix equations, is numer- 
ically quite intensive. 

Synthesis of non-symmetric PR gain: 
Algorithm 2 

For robust stabilization of passive systems, G does 
not have to be symmetric, and it suffices to have 
G + G T > 0. This condition is less restrictive than 
requiring symmetry, and therefore would generally re- 
sult in a smaller optimal value of the performance 
function. 


In [Moe.85], an algorithm was given for solving the 
standard LQ output feedback problem. It essentially 
involves solving two uncoupled Lyapunov equations 
at each iteration, and then restricting the step size in 
the direction of the resulting new value of G to en- 
sure closed-loop stability as well as reduction in J. 
The Lyapunov equations are linear, and therefore the 
computational requirement is quite reasonable. This 
algorithm can be modified for the design of PR con- 
trollers, by restricting the step size in the direction of 
new G (at each iteration) to ensure that G + G T > 0. 
This also ensures closed-loop stability at each itera- 
tion. Convergence cannot be guaranteed for this pro- 
cedure because the positivity constraint can drive the 
gain to the boundary. 

Numerical Examples 

Example 1: Flexible Structure Control 

The first example consists of a flexible space struc- 
ture with three lightly-damped elastic modes with 
frequencies (in rad/sec): 1.095, 2.3, and 2.6. The 
corresponding open-loop damping ratios are: 0.0023, 
0.0011, and 0.0019, respectively. Two actuators with 
collocated rate sensors are assumed, resulting in a 
weakly strictly positive real (WSPR) plant. An 
LQ performance index given by Eq. (6) is mini- 
mized by assuming the covariance of the initial state 
X 0 — /, and choosing the design variables Q = 
diag [10, 100, 10, 10, 10, 100], and A = hx 2 - 
Results with Algorithm 1: 

Algorithm 1 was first used for designing symmetric 
constant-gain optimal PR controller. The initial gain 
matrix was chosen to be diag[10, 10]. The initial value 
of J was 6.1398 x 10 3 . The minimization of J yielded 
an optimum value of J as 2.5276 x 10 3 . The resulting 
optimal symmetric, positive gain matrix was 

T 37.5665 59.9666 1 

~ |_ 59.9666 2 14.79 10 J ’ ^ 

The closed-loop eigenvalues of the system are given in 
Table 1. The position response at both sensor loca- 
tions are given in Figs. 1 and 2. It can be seen that 
the closed-loop response dies down within 30 seconds 
in both cases. 

Results with Algorithm 2: 

Algorithm 2, which permits non-symmetric gain 



matrix G , was next applied to the same problem, 
starting with the same initial value of G. The algo- 
rithm converged in 28 iterations, and the final value 
of J was: 2.5124e+03, virtually same as that obtained 
by Algorithm 1. The final optimal value of G was: 

j* 31.7064 36.7063 ' 

“ [64.3397 215.8499. 

The closed-loop eigenvalues, as well as the initial con- 
dition responses, were very close to those obtained by 
Algorithm 1. When different starting values of G were 
used, both algorithms converged to nearly the same fi- 
nal values. 

The example shows that both algorithms can ef- 
fectively design an optimal constant-gain positive- real 
controller. 

Example 2: Longitudinal Control System for 
F-18 Fighter Aircraft 

In the second example, linearized longitudinal 
models of an F-18 High- Alpha Research Vehicle 
(HARV) at four different flight conditions are consid- 
ered. The objective is to design a pitch- axis control 
system at 15,000 ft. altitude, and at the following 
combinations of speed and normal acceleration: (1) 
0.7 Mach and lg, (2) 0.6 Mach and lg, (3) 0.49 Mach 
and lg, and (4) 0.3 Mach and 0.37g. The control in- 
put is elevator deflection and the output is the pitch 
rate. This is the same system that was considered 
in [Kel.97], wherein the plant was first passified using 
a third-order series compensator with poles at —10, 
-.05, and —.0035 and zeros at —1, —0.5, and —0.08. 
This compensator robustly passifies the plant at all 
four flight conditions. 

The presence of the phugoid mode and the cor- 
responding zero causes numerical problems in control 
design algorithms. Therefore, passified short-period 
approximations (of order 5) were used for the purpose 
of controller design. The LQ performance function to 
be minimized has Q - / 5x5 and R — 10. The model 
corresponding to the second flight condition was used 
as the nominal model for controller design. 

Algorithm 1 was first used for designing a constant- 
gain PR controller for minimizing the LQ performance 
index J given by Eq. (6). The initial gain was chosen 
to be 0.1. The initial value of J was 1.09 x 10 4 . The 
minimization of J yielded an optimum with value of J 


as 483.3. The resulting gain was found to be 5.3264. 
The closed-loop eigenvalues for the nominal plant are 
given in Table 1 . The pitch rate responses for all flight 
conditions are given in Fig. 3. 

Algorithm 2 gave essentially identical results, as 
is expected because the controller gain is a scalar. 
The final compensator consists of the passifier in se- 
ries with the plant and the optimal gain in the feed- 
back path. The responses obtained were notably bet- 
ter than those obtained in [Kel.97] using LQG-optimal 
dynamic PR controllers of [Loz.90]. This could be at- 
tributed to the fact that only limited freedom is avail- 
able in choosing the performance function weights in 
the latter case. 

Gain Scheduling: The next step was to design op- 
timal controllers tuned to individual flight conditions. 
In each case, both algorithms converged to essentially 
same optima. The optimal gains for the four flight 
conditions are given in Table 2. The responses using 
individual optimized gains showed a slight improve- 
ment over those obtained using the optimal gain for 
the nominal design model (flight condition 2). If de- 
sired, the gains showed in Table 2 can be used for gain 
scheduling. 

Conclusions 

Synthesis of constant-gain positive-real LQ- 
optimal controllers was investigated for passive LTI 
systems. The controller design technique was demon- 
strated by two numerical examples. The synthesis 
methods presented in this paper, along with the ro- 
bust passification methods proposed in previous pub- 
lications by the authors, offer an effective tool for de- 
signing robust controllers for non-passive systems as 
well. 
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Table 1: Closed-loop eigenvalues 


Closed-loop Eigenvalues 

Algorithm 1 

Algorithm 2 

-0.1585 + 1.12181 
-0.6608 4- 2.17871 
-0.1627 + 2.5344i 

± 0.1370 + 1.1170i 
± 0.7000 + 2.1974i 
± 0.1587 + 2.5190i 
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Figure 1: Open-loop (dashed) and closed-loop (solid) 
response at sensor 1 



Figure 2: Open-loop (dashed) and closed-loop (solid) 
response at sensor 2 



Figure 3: Pitch rate response 


Table 2: Optimal gains 


Flight Cond. 

1 

2 

3 

4 

Optimal gain (G) 

5.56 

5.33 

4.33 

4.82 






